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1. Introduction 

The study of dissipation and decoherence is a central problem in the description of solid 
state systems both from the fundamental and the applicative point of view. One of 
the most popular models, employed to describe dissipation, is the one introduced by 
Caldeira and Leggett [lj in which the dissipative environment is a bath of harmonic 
oscillators whose fluctuations obey Gaussian statistics, linearly coupled to the system 
under consideration [21 13]. While this model is reasonable to describe the dissipation 
induced by large environments, such as equilibrium Fermi reservoirs, it is becoming 
clear that possible sources of decoherence do not necessarily fall into this category for 
nanodevices such as quantum dots and qubits. Solid state nanocircuits typically suffer 
from discrete noise, often originated from background charge fluctuations. Single/few 
impurities behaving either like random telegraph noise sources [4] or being entangled 
with the device have been recently observed (5] in different setups. The resulting 
non-gaussian effects can only be predicted introducing non-linear bath models |6j. In 
addition, the excellent control recently achieved in circuit-QED experiments [7] has 
paved the way to the observation of quantum effects originating from the coupling of 
multilevel nanodevices to a single harmonic mode of a high-Q cavity resulting in a 
effective bath having a structured power spectrum [Sj. All these mechanisms deviate 
from the canonical Caldeira and Leggett description and are examples of unconventional 
dissipative baths. 

The coupling to an unconventional dissipative bath arises also in the case of a 
nanoresonator coupled to a quantum dot, where electrons are coupled to few or even 
a single phonon mode [5]. An experimental realization of such a system is the phonon 
cavity studied by Weig et al. [10]. Here, the coupling to the single bath mode induces 
sizable effects on the transport properties of the dot, with the appearance of peaks in 
the conductance spectrum when the energy of tunneling electrons matches the energy of 
the phonon mode. A similar phenomenon occurs in a oscillating Ceo molecule contacted 
by gold leads [TTj . 

Often, in these cases both the leads bath and the localized bath mode are assumed to 
be in their thermal equilibrium [12]. However, there are several situations in which this 
assumption is not justified especially for the localized mode. For example, in suspended 
carbon nanotubes, the electrons are coupled with a single longitudinal stretching [13], 
or to a radial breathing |H] phonon mode. The unconventional single mode bath is 
driven out of equilibrium by the electron current and it induces the appearance of 
peculiar negative differential conductance traces [TSJ [IB] . Another example is that of 
the superconducting single electron transistor coupled to a single-mode nanoresonator 
studied by Naik et al. [T7] , where the back-action exerted by the single electron transistor 
cools down the phonon mode below the temperature of the environment. Then, it may 
be necessary to treat the leads and the single mode on a different footing. While the 
leads degrees of freedom can be assumed in their thermal equilibrium and are traced 
out in the usual way, the coupled dynamics of the electrons and the localized boson is 
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treated without assuming a priori a distribution for the localized mode. 

Such a issue will be the subject of investigation in this work. We will consider a 
spin degenerate single level quantum dot with finite Coulomb repulsion coupled to two 
external environments: 

(z) a pair of Fermi leads, described by a standard Caldeira Leggett model, 
(ii) an unconventional non-relaxed single phonon bath. 

The transport properties of such a system have been theoretically analyzed both in the 
sequential tunneling and in the cotunneling regime (T5J HEl EE] EUJ El]. The behavior 
strongly differs with respect to the one obtained for a system coupled to an equilibrated 
localized phonon mode [TB]. In the sequential tunneling regime, besides the negative 
differential conductance mentioned above, the zero-frequency shot noise S has been 
found to be even more sensitive to the bath properties [TH]. In the case of strong 
coupling to an out of equilibrium phonon bath, a super-Poissonian shot noise S > 2el 
(— e is the electron charge, I the corresponding current) has been predicted [19]. On 
the other hand, for a thermal bath, sub-Poissonian shot noise S < 2el is always 
found [20]. In the cotunneling regime, vibrational absorption sidebands occur within 
the Coulomb-blockaded regime, which disappear for the case of a single bath in thermal 
equilibrium [T9l |2T| Q3] • This further confirms the importance of such a unequilibrated 
phonon bath to describe experimental systems. 

In order to consider the coherent dynamics of the quantum dot with the single 
mode bath one can derive a generalized master equation within the Born-Markov 
approximation for the reduced density matrix of the system and the localized bath, 
by tracing out the leads bath. Such an approach has been employed already in the 
past for describing systems such as a metallic [22] or a superconducting single electron 
transistor [23] . 

Since the transport properties have been already studied in great details, we will 
concentrate on the steady state out of equilibrium phonon distribution of the bath 
(ii) induced by the tunneling electrons, which in the past has received much less 
attention [TSJ [2U [22J EZ]- In particular, we will analyze the phonon Fano factor 
var(/) _ (Z 2 ) - (Z) 2 
ph (I) (I) ' 1 ' 

defined as the ratio between the variance of the phonon occupation number var(Z) and 
the average occupation number (I). The Fano factor brings information about the 
statistics of the single phonon bath mode. For a thermally equilibrated bath, one always 
obtains a super-Poissonian Fano factor F P h > 1. Such a super-Poissonian value is typical 
for "classical" boson distributions, as for example for photons in classical light [28J. 
As it will be shown later, in most of the transport regimes, out of equilibrium bath 
distributions arise which, although being strongly different from a thermal distribution, 
display a super-Poissonian character. However, it is possible to find particular transport 
regimes in which F p h < 1, corresponding to peculiar bath distributions [23l|27]. This case 
is analogous to that of the reduced photon fluctuation found in quantum non-classical 
radiation [291130]. 
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Our main results are the following. In the regime where the single bath mode 
oscillations are much faster than the average electron dwell time, we adopt the rotating 
wave approximation and we investigate the phonon Fano factor. We identify the 
transport regimes where a sub-Poissonian single mode bath occurs for the case of a 
single occupancy of the dot. For a double occupancy, additional transport regions in 
which F p h < 1 occurs. We evaluate the "stability boundary" of the sub-Poissonian 
bath as a function of tunnel barriers asymmetry and coupling strength with the bath 
mode, finding that F p h < 1 is destroyed for increasing damping of the single mode bath. 
Finally, we present some preliminary results on the effects of the coherent phonon bath 
dynamics occurring when the oscillation time of the localized mode is of the same order 
or smaller than the electron tunneling time. 



2. Model and methods 



2.1. The system and its environments 



We model the system as a small quantum dot, with an average level spacing of the same 
order than the charging energy, described as a spin-degenerate single level with on-site 
Coulomb repulsion [31] 



H s = ^n- en^V g + |n(n - 1) . (2) 

Here, £ is the single energy level, n = Yla=± n <* * s ^ ne total occupation number with 
n a = d\d a the occupation of spin o = ±1 (units h/2), and d a {d) a ) are the fermionic 
annihilation (creation) operators. The second term in Eq. (|2| is the coupling to an 
external gate voltage V g , with — e the electron charge, C g the gate capacitance, and 
C = e 2 /U the dot capacitance. Choosing as a reference level £ = U/2 one has 
H s = en + Un(n — l)/2 with e = U(l — 2n g )/2 written in terms of the number of 
charges induced by the gate n g = C g V g /e. Note that the choice of £ determines the 
values of n g for which resonance between the n, n + 1 states occur: a different choice of 
reference for £ would simply result in a shift of n g . 

The system is coupled to two dissipative baths with Hamiltonians (h = 1) 

H^=u Q tfb- iff = ?k cl Aa c aAa . (3) 



a=L,R,k,a=± 



Here, H^ 1 ' describes the harmonic oscillator single bath mode with frequency ujq, and 

(2) 

HI represents the external left (L) and right (it) leads of noninteracting electrons 
with fermionic operators c a ,k,a and v aka . The localized phonon mode is undamped, a 



discussion of possible damping effects is deferred to Sec. |4.1| The leads are assumed in 
thermal equilibrium with respect to their electrochemical potential Hl,r — A*q ± eV/2, 
where V is the applied voltage, and fj,Q is the reference chemical potential. The 
noninteracting Fermi leads can be mapped onto a ohmic dissipative bath within the 
Caldeira Leggett formalism [3j. More detailed models, involving a anharmonic (23] or 
distorted [25] oscillator or the presence of interacting Luttinger liquid leads J2H], which 
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have been studied in the past especially concerning their transport properties, will not 
be addressed here. 

The coupling of the system with the single-mode bath is linear in the oscillator 
coordinate and in the effective charge n — n g on the dot 

= V\u (b + tf)(n-n g ). (4) 

This describes realistic situations such as a single electron transistor capacitively coupled 
to a vibrating, charged gate [321 EH] • The tunneling Hamiltonian 

H T = E *« E c L«a + h - c - • ( 5 ) 

represents the coupling of the dot with the external leads, through tunneling amplitudes 
tz and £r. In the following, a trace over the leads degrees of freedom will be performed. 

The full dynamics of the single mode will be retained, and its degrees of freedom will 
not be traced away. It is possible to diagonalize exactly the Hamiltonian representing 
the dot coupled to single phonon mode 

H X = H S + + H$ , (6) 

by means of the Lang-Firsov polaron transformation [IB], with generator U = 
exp [— rj{b — fat)], and r\ = V\(n — n g ). The transformed operators, denoted by an 



while n = n. Upon 



overbar O = UOU.\ are b = b — rj and d a = d a exp \/~X(b — tf) 
transformation, we have 

H x = + ^n(n - 1) + u l , (7) 
with renormalized level position e and Coulomb repulsion U 

e=|(l-2n fl ); U = U - 2Xu . (8) 

The eigenvectors of Eq. (J7| are denoted by \n,l), with energy E n j. Note that the 
polaronic renormalization of the Coulomb interaction may lead to two qualitatively 
different physical scenarios for U ^ or equivalently A ^ U/2uq. The regime U < 0, 
where the single occupation is always forbidden and the sequential transport is blocked, 
has been considered by several authors [3U [351 [36] . I n the present paper, we will consider 
the opposite case with U > treating the sequential tunneling regime. 
Finally, the transformed tunneling Hamiltonian is 

= E ^ (6 ~ 6t) E + h - c - > ( 9 ) 

a=L,R k,a=± 

with an the explicit dependence on the oscillator variables. 



2.2. Generalized master equation 

The dynamics of the dot and single bath mode can be described in terms of the reduced 
density matrix p(t), defined as the trace over the leads bath of the total density matrix 
Ptot{t)'- p{t) = Tii{p tot (t)} . In the following we will mainly work in the polaron frame 
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with the density matrix denoted by pit)-, the evolution in the original frame being easily 
traced back via the canonical transformation p(t) = W(t)p(t)U(t). 

We will consider the regime of weak tunneling, treating Sj^ in Eq. (ji) to lowest 
order, and assume the characteristic memory time of the leads much shorter than the 
response time of the dot interacting with the localized phonon bath. This allows to 
treat the dynamics in the Born-Markov approximation [37] leading to a generalized 
master equation (GME). Assuming a factorized total density p to tif) = p(t) <S> Pi(0) with 
the leads in thermal equilibrium with respect to corresponding chemical potential Pl,r- 
Pl(0) = Pl(0) <S> Pr(0), the time evolution in the interaction picture (denoted with the 
subscript "/") is 



P~i(t) 



J2 / drUQjA^QUt-^p^K^T) 



[QiMMVQPiA 1 ~ r )^» + [QlMQiA* - T)pj(t)]K-(-r) 

- [Q{ a (t),Pi(t)QiAt-r)}K + (-r)} . (10) 

Here, we defined the operator Q a = e^ 6-6 and the leads correlation functions 

K ± (r) = J2\t a \ 2 e iekT f^(e k ), (11) 

a.k 

with f£(e) = 1/{1 + exp[/3(e — p a )}} the Fermi distribution of lead a and f~(e) = 
1 — fa{e). Note that the correlation function K + (t) can be cast into the form [38j |39j 

K+(r) = iyJ2 \t a \ 2 e l ^ T uJ c e- wiT) , (12) 
a 

where v is the leads density of states and 

W(t) = J™ du ^ | [1 - cos (ut)] coth f^j + i sin(cur) J , (13) 

represents the dissipative kernel of a set of harmonic oscillators with ohmic spectral 
density J(u) = coe~ UJ ^ ulc , and frequency cutoff u c [3]. Expression (12) clarifies the 
connection between the electronic leads and the bosonic modes of the Caldeira Leggett 
model. 

The above approximations allow to describe the dynamics of the dot in the so called 
sequential tunneling regime. They are typically valid for temperatures larger than the 
level broadening induced by tunneling /c#T ^> I\ + T^, with r Q = 2itv\t a \ 2 . At lower 
temperatures, coherences between the leads and the dot play a crucial role [18J. These 
effects will be not discussed in this work. 

It is convenient to project Eq. (10) on the eingenstates of the Hamiltonian ^ and 



perform the time integrations with the aid of J °° dr e lEr = ir5(E) + iP.V.(l/E). In the 
following, terms stemming from the principal value are neglected since they lead to small 
corrections in the perturbative regime [221 HQ]. Denoting with p™ qq ,(t) = (n, q\pj(t)\n, q') 
the GME has the compact form (here and in the following, < n < 2 unless stated 
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p,p 



0<n+fe<2 

~~ z ^ A p',p 7pp> [ A p 'g PiM'tfje + A p>qi P/, gP W e qp \(- l i4 J 



Here we defined uj pp i = uj ■ (p — p'), with the indexes p,p',q,q' that run over [0, oo). 
The factors z = z 2 = 2, z\ = 1 arise due to the spin degeneracy of the state n = 1. 
Furthermore, it is 
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,n,n+l _ f + / p __ p 



r 



7, 



v n,n+l v 



and we have introduced 

with = (n, q\e^( b ~ b ^ \n, q') the Franck-Condon factors, given by 

W-q\ 



X qq i — e 



X, \ —A'-%\) \ Sg a( q '-q)VX 



(16) 



(17) 



where q < = min{g, q'}, g> = m&x{q,q'} and L V (X) is the generalized Laguerre 
polynomial. They satisfy the relation X% , = . 



In the Schrodinger representation, Eq. rtl4h becomes 



lU q qit 



Tin (+\ i \ ^ \ ^ J „ v n+k,n v n+k,n/ „n+k,n , „ n+fc,n\ -n+fc /+\ 
Pqq'i t ) + 2^ l^\ Z ^+k X q'p' X qp ilpq + %>q> ) Ppp' W 



T^n,n+fc n,n+fc r y n, n+fc -n /,\ . v n,n+k -n \ 
Z « A p'p 7 pp ' [ A p'q Ppq'^) + A p>q> Pqp^llj 



The GME (18) is an infinite set of coupled linear differential equations which cannot 



be solved analytically under general conditions. We numerically solve the system 
truncating the harmonic oscillator Hilbert space to increasingly larger sizes q, q' < 70, 
until convergence is achieved. 



2.3. Rotating wave approximation 

In the regime of fast vibrational motion of the localized bath, uq ^> r^, Tr, one can 
perform the rotating wave approximation (RWA), neglecting terms in Eq. (14) with an 
explicit oscillatory exponential time dependence [37]. In the Schrodinger representation, 
this leads to a GME of the form 

( oo 

—n i ,\ _ • -n , \ ^ J \ ^ x^n+fc,n -yn+k,n -n+k / n+k,n . n+k,n \ 

Pqq'K 1 ) — ~ lU, qq'Pqq'+ /-^l \ 2-^1 Zn+k 1 9,?+Z> q' ,q'+pPq+P,q'+P \lq+P,q +7 g '+p,q'J 



k=±l {p= — q < 

0<n+k<2 



-Znp n 



oo 

E\ ( v n,n+k\2„n,n+k . / v 
[( X pq' ) Iq'p + ( X pq 



n,n+k\2 n,n+k 
I Iqp 



p=0 



(19) 
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Note that in the RWA the density matrix elements p™ , and p$ are coupled only if 
V — I = q' — q. This implies that the diagonal elements pq q {t) are decoupled from 
the off-diagonal ones. In addition, the non diagonal elements vanish in the stationary 
regime, p" z ,_^(t — > oo) — > 0, as a difference with the coherent regime (see section 4.2). 
Hence, within the RWA all the stationary properties are determined by the diagonal 



occupation probabilities p qq = P nq and Eq. (19), in the stationary limit, assumes the 
form of a standard rate equation 

oo oo 

z n p nq e EC- E £ww£-* ,B = o ( 2 °) 

fc=±l p =0 k=±i p=0 

0<n+fc<2 0<n+fc<2 

pnn' -pnn' , -pnn' cy -yl n J ln ' (0~\ \ 

pq L,pq i" R,pq p,q Ipq \ ) 

with Tp™' the tunneling rates for the transition \n,p) — > \n',q) and the tunneling 
rate on barrier a. We specify here the energy regions where tunneling processes with 
rates T™' are allowed (i.e. T™' ^ 0) for temperatures low enough that the Fermi 



function can be approximated as a step, T u^/ks- From the definition (15) follows 



Kq +1 + => E n+lA ~ E n,P ± < (22) 



pq 

r? 1,n ^0 =► E n+ljP -E njq ±eV/2>0, (23) 
with the + (— ) sign for a tunneling event on the right (left) barrier. The Franck Condon 



factors in the tunneling rates Eq.(21) are responsible for two important effects: 

(i) suppression of all the tunneling rates by a factor e~ A , 

(ii) non-trivial dependence on the phonon indexes p, q. 

In particular, for moderate interactions A < 1, transitions which conserve or change 
slightly the phonon number have the largest rate and, among these, the ones involving 
small p and q are dominant. Vice versa, for A ^> 1, transitions with a large change 
in the phonon number are favored. These considerations will play an important role 
in explaining both the transport properties and the characteristics of the single mode 
bath. 

Finally, the stationary current can be expressed in terms of the occupation 
probabilities as follows 

I = ~e E E M^Iw 1 - z n+1 P n+1>l rZ¥] , (24) 

n=0,l q,q' 

note that I is independent on the position. 

In order to obtain the stationary solution of the phononic populations, we numerically 



solved Eq. (20 ) for a sufficiently large set of harmonic oscillator states with q < 60, which 



guarantees a fair convergence. A discussion of the range of validity of this approximation 



will be presented in Section 4.2 



3. Results 



In the following we will present several effects induced by the coupling of the quantum 
dot with the unequilib rated single bath mode. We will consider the regime /c#T u 
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Figure 1. (a) Current / (units — cTl) as a function of V (units uoq/c) for an 
unequilibrated (solid) or a thermal (dashed) phonon bath, for n g = 0.525, A — 0.01, 
A = 3, k B T = 0.01 uj , U = 20 w and T L = 10~ 3 w . (b) Phonon distributions P h 
as a function of I, for the cases of a thermal bath (blue) and an unequilibrated bath 
calculated for eV — 2 uio (red) and eV — 8 uj (green), all other parameters as in panel 
(a). 



the current /, numerically evaluated from Eq.(24) in the RWA, as a function of voltage in 
the presence of an asymmetry A = |£r| 2 /|£z.| 2 < 1. Two different conditions of the single 
bath mode are considered: a completely out of equilibrium distribution Pi, determined 



by solving the rate equations ( 20 ) , and a thermal distribution P^ 



P, 



(th) 



-lf3u>o , 



1-e 



-f3uo 



(25) 



and j3 = (/cbT) -1 . For V < u /e, the dot is in the Coulomb blockade regime. For 
increasing V, steps in the current reflect phonon excitations in the bath. The solid line 
is the result with Pi, the dashed line displays the current obtained imposing a thermal 
distribution of the bath mode. The corresponding phonon distributions Pi deviates 
considerably from -P/*^, as shown in Fig. l](b). While P^ at T uJo/ks is almost 
concentrated at / = 0, the unequilibrated distributions are broader. 

The behavior of the current is qualitatively different, depending on the bath 
conditions. In particular, for V < 5 u^o/e the current in the equilibrated case is 
smaller than for an unequilibrated phonon bath while for V > 7uo/e the situation 
is reversed [18, 20J. This fact can be explained as follows. For A > 1, as in the case 
of Fig. [IJa), the dominant transition involves / = — ► int(A) [20j, as a consequence 
of the Franck Condon factors (17). For the parameters in the figure, it corresponds 
to the transition I = — > 3, which is open only for V > 5co>o as can be checked by 
inspecting Eqs. (22), (23). A large occupation of the states with I = leads therefore 



to an increase of the current. This explains why, for large voltages V > 7ooo/e, in 
the case of a thermal distribution P^ one obtains a larger current with respect to 
the unequilibrated case where a broader phonon distribution implies a smaller P [20J . 
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Vice versa, for V < 5 u>o/e only small tunneling rates with Al < 2 are open, therefore 
the larger the number of transport channels, the larger is the current. In the case of 
a thermal bath only transitions originating from I = can contribute to the current, 
while in the unequilibrated case also transitions starting from / = 1 and / = 2 states 
contribute, resulting in a higher current with respect to the equilibrated case. 

Also the current fluctuations display deviations (not shown here): in the 
unequilibrated case a super-Poissonian current noise appears, while for a thermal 
bath one obtains a sub-Poissonian noise [T9J 120] . These topics have been studied in 
great details [151 HH1 EE |2T], thus we will not repeat here the analysis of transport 
properties induced by the unconventional localized mode. 

Instead, we will focus on the properties of the single bath mode, and we will characterize 
the out of equilibrium phonon distribution induced by the current flow. In particular, 
we will analyze the Fano factor. The averages involved in Eq. (fl]) are equivalently 
defined in the original frame or in the polaron one as (O) = Tr{Op} = Tr{(5p} with the 
trace performed over the degrees of freedom of the system and of the single bath mode. 
Since operators in the polaron frame are expressed in terms of the ones in the original 
frame, it will be particularly useful to introduce the hybrid average (0) p = Tr{Op}. 
With respect to this average the above occupation number and variance are given by 
[rj = y/x(n - rig)) 

(I) = (l) p - ( V x) p + ( V 2 ) p (26) 
(l 2 ) = (1 2 )p + (*% + faV), + 2{rfl) p - 2{rfx) p - ( V {1, x})- p , (27) 

where {a, b} = ab + ba and x = b + w. Note that the Fano factor can also be connected 
to the Mandel parameter Q = F ph — 1 [JT] , introduced in quantum optics to discriminate 
the sub-Poissonian statistics and photon antibunching [291 ED] • 

We point out here that in the case of a thermally equilibrated bath, with a diagonal 
density matrix Eq. (25) right side, one has F p h > 1 [|J As we will see below, deviations 



from an equilibrated bath usually result in super-Poissonian distributions. However, 
in special transport regimes it is possible to obtain peculiar distributions with a sub- 
Poissonian character. 

3.1. Single occupancy 

We start from the case U ^> u within the RWA. In this regime the dot, around 
the resonance condition n g ~ 1/2, has a single occupancy with n = 0, 1 and shows 
qualitatively similar behavior to that discussed in Ref. [27], where U — > oo. Within the 



RWA, the Fano factor is simplified since some averages in (26) and (27) vanish due to 
the diagonal form of the density matrix. It is convenient to express F p h in terms of the 
numerator of the Mandel parameter F ph = 1 + Q num / (I) where 

Q num = Q- P + Mv 2 1)p ~ 2(v 2 )p(1)p + (V 4 )p ~ (V 2 )l (28) 



I Note that in the thermal regime, also the density matrix in the original frame is diagonal if expressed 
on the basis of the eigenstates of H\, given by W\n, I). 
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Qp =(l 2 h-(l)l-(l)r- (29) 
Here, we have separated a contribution, Qp, depending only on the phonon distribution, 
and terms related to charge fluctuations on the dot (terms in if and r] A ). In general, 
one has Q n um ^ 

and Q num > Qp |27j . The necessary condition to have F v h < 1 is 
Qp < 0. However, this does not guarantee sub-Poissonian behavior since the variance 
(r] 4 )p — (r] 2 )p, being definite positive, could drive F ph > 1. 

The quantity Qp in (29) is expressed in terms of the phonon distribution Pi in Eq. (25) 
as follows 

Q,=x^z(i-i)-feflA • ( 3 °) 

A sub-Poissonian distribution can originate from peculiar phonon populations 
characterized by Qp < 0. This may be achieved under the following selective condition 

P ~Pi>iW (31) 

Phonon distributions satisfying the above condition in the following will be referred 
to as "selective populations". They will play an important role in inducing a "non- 
classical" phonon gas. We now analyze the parameter regime where this condition can 
be achieved. 




Figure 2. (a) Color map of the numerically calculated phonon Fano factor K,/, as a 
function of V (units coo/e) and n g for A = 0.1, A = 0.8, ksT = 0.01 u>o, U — 20 loq and 
= 10~ 3 wo- (b) Main: F p h as a function of V (units u>o/e) calculated along the line 
U(l — 2n g ) + eV — ojq — in the V, n g plane (corresponding to the black arrow in the 
bottom right scheme). Inset: phonon distributions Pi as a function of I for different 
values of A, calculated in the middle of region Ri (n g — 0.525, eV = 3ojq); other 
parameters as in (a). For A = 0.1,0.8,3 (left to right panels) it is F p h « 1.2,0.9,1.3 
respectively. 



Fig. |2^a) shows a color map of the numerically calculated F p h as a function of V 
and n g for U = 20 u . It displays a checkerboard pattern where each region is 
characterized by the activation of transport channels involving transitions between 
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different phonon states. The regions R p and R* will be particularly important for the 
following discussions. Their center is at 

= ; n<f = l±^ P (32) 



with the + (— ) sign for region R p (Rp). As stems from (22), (23) the following conditions 
are necessary for transitions to be open 

\0,l + q) -»• + 2u (l' -l) + U(l-2n g )±eV <0 (33) 

\l,l + q) -»• \0,l' + q) & 2co (l' -I) -U(l-2n g )±eV <0. (34) 



For the values of V and n g within region R p , it follows that the transitions (33 ) are open 
if I' < I + p, while for the transitions (34) one needs V <l. Similarly, in regions R* one 
finds that for the transitions (33) one has /' < I while for (34) it is I' < I +p. 
Thus, in region R the transition rates towards excited phonon states are all closed, 
therefore only the I = phonon state is populated. Here, the solution of the rate 
equations is simply given by P o = A/(A + 2), P 10 = 1 — P o (f° r T -C ujQ/k B ) and one 
gets a super-Poissonian behavior 

2A (2n g - l) 2 

Fph ~ 1 + X 2 + A(2 + A)nj-4n g + 2 > 1 ' 

The situation is more interesting in the other regions displayed in the figure. Indeed, in 
most of them the out of equilibrium phonon bath displays a super-Poissonian character. 
However, one can note that sub-Poissonian values, F ph < 1, are present at low voltages 
in region R\ as shown also in the main panel of Fig. [2](b) . Concentrating on Ri, the 
selective population arises if the following conditions are met: 

(i) the allowed phonon transitions \n, I) — > \n', I') must satisfy V < 1 + I, which implies 
no direct population of the states with I > 2 from an initial state |n, 0); 

(ii) the asymmetry has to be A < 1, which suppresses the rate for tunnel-out transition 
|1, 1) — > \0,l') since for V > electrons flow from the left to the right lead. 
Condition (z) is precisely realized in region R±. By virtue of this, in order to populate 
excited phonon states with I > 2 starting from the phonon ground state I = within 
Ri, one needs to perform at least the following sequence of transitions 

|0,0)-|1,1)-|0,1)-|1,2)-... • (35) 

This implies that at least I— 1 tunnel-out events are performed in R\ in order to populate 
an excited state with / > 1 and this leads to Pi>\ oc A 1 ^ 1 . 

Condition (ii) and the above result induce a small occupation of the / > 2 states for 
A < 1 [27J. The inset of Fig. |2^b) shows the phonon distribution in region Ri for 
A = 0.1 and three different coupling strengths. One can observe that all distributions 
are markedly out of equilibrium and very different from a thermal population, see the 
blue bar in Fig. [T|b). In addition, only the case of intermediate coupling A = 0.8 
(second panel) displays a selective population as defined above. For small A (first panel) 
broad phonon distributions are obtained. On the other hand, for A > 1, the phonon 
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distribution gets very narrow with Pq ^> Pi>\ and no selective phonon population as in 



(31) may be achieved. 

As a result of the competition between asymmetry, coupling strength and charge 
fluctuations, it is possible to find F p h < 1 only in the parameter space characterized 
by both A < 1 and A < A c (\) < 1 where A C (X) is a threshold asymmetry of order unity 
for A ~ 1 [27] . Furthermore, it is possible to show [27] that similar arguments apply in 
the region S with a much stronger asymmetry, while region R*[ is similar to Ri, inverting 
the asymmetry: A — > 1/A. 

In all other regions R p > 2 and R* p>2) the transitions \n, I) — > \n',l') with I' < I + p 
and p>2 are active. Now, the transitions \n, 0) — > \n',l') with V < p directly populate 
phonon states with V > 2 from the phonon ground state \n,0). As a result, the 
mechanism to obtain a selective population discussed above is no longer effective, and 
broad phonon distributions are attained for A < 1. They neither display a selective 
population, nor F p h < 1. Finally, for A > 1, narrow distributions arise and again no 
F p h < 1 is obtained. 

3.2. Double occupancy 

We now turn to analyze the effects of the dot double occupancy on the phonon 
distribution. We focus on the case of A < 1 and A < 1 which, as we discussed above, 
favors F p h < 1 when is large. The choice A < 1 is also favorable for obtaining 
double occupancy, as tunnel-out rates are suppressed. We will concentrate on the 
regime > 2ujq that is relevant in realistic situations with an intermediate dot-phonon 
coupling strength. Figs. [3|a)-(c) show F ph for the same parameters as in Fig. [2|a) but 
much smaller 0. In addition to region R±, it is now F p h < 1 in region R% for = 5 ujq 
and in region R 2 for = 4 uj q . A plot of F ph calculated along the diagonal line shown in 
Figs. [3|a) andj^b) is displayed in panel (c). Note that the new regions where F p h < 1 
shift at larger voltages with increasing U. This explains why the case U ^> u> is similar 
to the case U —>■ oo near the resonance n — — ► 1, where only the region Ri displays 
Fph < 1. The above behaviors represent special cases of the non-trivial general result: 
for obtaining F p h < 1 in region R p > 2 it is necessary to have 

(p + 1 Vo < U < (p + 3 Vo . (36) 



In the following, we will explain the origin of Eq. (36) and show that it relates to specific 
transitions which involve the state n — 2. 

Let us first of all determine the conditions under which double occupancy can be 
achieved. In the absence of phonons (A = 0) the onset of the transition n = 1 — >• 2 occurs 
for V > Vo = U(3 — 2n g )/e. On the other hand, in the presence of phonons, the state 
1 2, 0) can be occupied for smaller V: indeed, the transitions |1, l + q) — > \2,q) with I > 1 
are active for V > Vi = [0(3 — 2n g ) — 2lLU ]/e, with V\ < V . This phonon-mediated 
mechanism for obtaining a double occupancy is similar to that described by Shen et 
al. (12] for the occupation of the LUMO in a molecular quantum dot. The condition 
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Figure 3. (a-c) Color map of the phonon Fano factor F p h in the V (units oJo/e), n g 
plane for A = 0.1, A = 0.8, fc B T = 0.01 uj q , T l = 10~ 3 w and different values of U: 
(a) U — 5 wo, (b) 0" = 4 u , (c) f7 = 3 uj . (d) F p /, as a function of V (units cj /e) 
calculated along the line U(l — 2n g ) + eV — cj = for U — 5 ojq (solid) and U — 4 uj 
(dashed) (see the black arrows in panels (a), (b)). 



for having the transition |1, 1 + q) 
relevant for V > 0) is 



\2,q) open on the left barrier (forward transition, 



U(l 



2n g ) — 2lu 1 



2U - eV < . 



(37) 



In view of voltages constraints (32) for region R p , the above transitions open there for 
U < (p + 1 + l)o>o- We will now show that double occupancy can lead to a selective 
phonon population with P ~ Pi P/>2 in the regions P p >2- For U < (p + 3)uj the 
transitions 

\l,l + q) ->\2,q) l>2 (38) 

are active within R p and decrease the occupation probability of the excited phonon bath 
states with I > 2. They come in addition to the usual transitions 

|1,Z)-HO,0 l' <l (39) 

which also depopulate the excited phonon states and were already present when U ^> ujq. 
Note that these latter transitions alone are not sufficient to give rise to a selective phonon 
population: indeed, broad phonon distributions are obtained for A < 1 in R p >2 as already 



discussed in Sec. 3.1 Now, since for A < 1 and V > tunnel-in events are faster than 



tunnel-out ones, the new transitions (38), are much more effective than transitions (39) 



in depopulating states with I > 2. With these new fast channels it is now possible to 
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obtain lower occupations of all the excited phonon bath states up to / 
leading to a selective phonon population. 



2, eventually 



This explains the right hand side of the inequality (36). To obtain the left hand side 



of (36) it is sufficient to recognize that if U < (p + 1)uq the additional depopulation 
transitions |1, 1 + q) — > |2, q) with / > 1 are open in all R p . They also deplete the states 
with 1 = 1, eventually leading to too a narrow distribution with Pq « 1 and F p h > 1. 
Therefore, tuning U it is possible to trim the shape of phonon population from broad, 
to selective, up to a very narrow distribution. Figure [| shows Pi as a function of I, 




Figure 4. Phonon distribution Pi as a function of I calculated in the center of region 
i?2 for different values of U, all other parameters as in Fig.J^a). For U = 5ujo, Aujq, 3ojo 
it is Fph f» 1.1,0.8,1.5, respectively. 



calculated in the center of R2 for different U, all other parameters are the same as in 
Fig. [2^a). For U = 5u , the transitions \ l,2 + q) — > |2, q) are closed and a broad phonon 
distribution (blue bars) with F p h > 1 is obtained. When U = 4 u , the above transitions 
are open and the channel |1,2) — > |2,0) creates the selective population displayed by 
the green bars. Further reducing U = 3 Uq, also the transitions |1, 1) — > |2,0) opens up 
and the resulting narrow distribution displays super-Poissonian character. 



4. Damping effects and the validity of the RWA 



In Section 2.3 we have discussed the out of equilibrium localized bath arising in 
the transport regime, concentrating particularly on the occurrence of sub-Poissonian 
phonon distributions. The results were derived assuming that the localized phonon 
mode has an infinite life time in the absence of tunneling. In addition, up to now we 
have always considered the case of a fast motion of the oscillator in comparison to the 
transport time scales and adopted the RWA. In this section we study the robustness of 
the results, when the above conditions are relaxed. 
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4-1. Towards a relaxed single mode bath 



In order to damp the single mode bath towards an equilibrium distribution, we couple it 
to a conventional dissipative environment, represented by a set of harmonic oscillators 



H, 



(3) 



(40) 



with a linear coupling 



(41) 



and bath spectral density J{ui) = 2nJ2j UJ ]x'j^( u ' ~ ty). 111 the presence of the 
additional bath (40), it is still possible to diagonalize exactly H\ by means of a canonical 
transformation which also involves the operators bj, b^ [P2| [20] . Following a derivation 



analogous to that discussed in Sec. 2.2 and treating both Eq. ^ and Eq. (41) to 
the lowest perturbative order, we obtain the following stationary rate equations in the 
RWA [121 [20] 

00 00 

_ y p \ r \ r pn,n+fe , \ r \ r „ p pn+fc,n 
&U 1 nq / J / J 1 qp T / j / j ^n+k 1 n+k,p L pq 



fc=±l p=0 

0<n+fc<2 



fc=±l p=0 
0<n+fc<2 



p=0 

The relaxation rates 



n V^ pq 



p=0 



nq 1 - qp 



0. 



-prel 
1 9,9-1 



' 1 9-1,9 



qwT L [l 



y 



(42) 



(43) 



induce transitions \n, q) 



\n, q ± 1) and drive the single mode bath towards thermal 
equilibrium, while the tunneling rates T™' are those described in Sec. 2.3 Here, 
w = J(uq)/Tl parametrizes the relaxation strength. For w — > oo, phonons are 
completely relaxed and Pi — > Pf- , characterized by a super-Poissonian Fano factor. 
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Figure 5. (a) Phonon Fano factor F p h for A = 0.1 (solid) and A = 0.01 (dashed) as a 
function of the relaxation strength w (log scale) for [7 = 4 ujq, eV = 3 lu and n g = 0.75 
(the center of i?2). All other parameters as in Figs. |3ja)-(c). (b) Phonon populations 
Pq (solid), P\ (dashed) as a function of w (log scale) for the same parameters as in (a). 



The effects of relaxation on the phonon distribution are illustrated in region R 2 for 
= 4: ujq. All other regions display qualitatively similar behavior. Fig. [5Fa) shows F p h as 




0.01 0.1 1 A 10 0.01 0.1 1 A 10 



Figure 6. (a) Color map of F p h as a function of A, A (both in log scale), for eV — 2 u>q, 
n g = 0.625 (region Ri), U — 4 w and w — 0.01, other parameters as in Figs. [3ja)-(c). 
The black contour line is set at F p h — 1. (b) Same as panel (a) but with w = 0.1. 



a function of w, calculated in the center of i?2- Asto increases, F p h crosses from F p h < 1 
in the unequilibrated regime to a super-Poissonian maximum, before reaching the 
thermal value which, for T <C uJo/ks, is only slightly above 1. The effects of relaxation 
become relevant when the typical life time of the excited phonon bath T re \ = Jiuoo)^ 1 
is comparable with the slowest typical charge dwell time r tr = (min{l, A}YLe~ x )^ 1 . 
This explains also the shift of the maximum towards smaller w for smaller asymmetries 
observed in Fig.[5|a). In Fig.[5|b), one can see that for increasing relaxation, the selective 
phonon population is destroyed with a tendency towards an almost full occupation of 
the I = state, in accordance with a low temperature (see Eq. (25). 

We now analyze the robustness of the sub-Poissonian phonon bath against w 
varying the parameters A, X. Figures[6^a) and[6^b) show a color map of Fph as a function 
of A and A in the middle of the region Ri for two different values of the relaxation 
strength. Superimposed to the color map, the black contour line signals F ph = 1. This 
line corresponds to the "stability boundary" of the sub-Poissonian phonon distribution 
with respect to asymmetry and coupling strength. Increasing relaxation, Fig. ^b), the 
region of A and A where F p h < 1 shrinks, the regime A <C 1 is more strongly affected, 
in accordance with the discussion above, while the sub-Poissonian distribution is still 
stable for not too small A < 1 and A ~ 1. Qualitatively analogous results are observed 
in every region R p . 

4.2. Validity of the RWA 



When the dynamics of the single mode bath is not much faster than the typical electron 
tunneling rate the RWA is no longer justified and the fully coherent GME in Eq. (18) 
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must be solved. In this case, also the non diagonal elements of the density matrix 
will be different from zero in the stationary regime. Even more important, while in the 
RWA the diagonal elements p% of the density matrix are decoupled from the off-diagonal 
ones, in general this is no longer the case. Thus the coherences pip^i may influence the 
phonon distribution P n i. In this last part, we will present preliminary results showing 
this effect. For simplicity, we consider U ^> uq, without relaxation w = 0. 
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Figure 7. (a) Solid: phonon Fano factor F p h obtained from the solution of Eq. (18 1 
as a function of 7 = T l /uj for A = 0.01, A = 0.1, T = 0.1 uj /k B , U = 20 lu , and 
eV = 2 ujq with n g = 0.525 (red, in region R\) or eV = 3 ujq with n g — 0.5 (blue, in 
region S). Dashed: F p h obtained from the RWA. (b) Relative difference A between F p h 
and F® h (see text for definition) , as a function of 7 for the same parameters as panel 



(a), (c) Diagonal elements Pi of the solution of Eq. 18 as a function of I calculated for 



eV = 2 uj and n g — 0.525 and different values of 7. (d) Same as (c) but for eV = 3 uj 



and n g = 0.5. 



Fig. shows F p h obtained from the solutions of the full GME in Eq. ( 18 ) as a function 
of 7 = Tl/ujq, where we keep u fixed (note that the calculation of F ph via Eqs. (26) 
and (27) contains also averages of non-diagonal operators). The red curve corresponds 
to Fph in the center of region Ri, while the blue one gives F p h in region S. The dashed 
lines correspond to the RWA solution. In both cases, it is clear that deviations from the 
RWA (dashed line) occur for 7 > 0.05. For increasing values of 7, F p h saturates into a 
new regime, characterized by sharply reduced phonon fluctuations. 

In order to assess the relevance of the off-diagonal density matrix elements, we have 
evaluated the ratio A = \F ph — F® h \/F ph , where F® h is calculated from the GME solution 



Phonon distributions of a single bath mode coupled to a quantum dot 



19 



neglecting the coherences Pii^i- We show A in Fig. |7](b): clearly, the impact of the 
off-diagonal terms on F p h is very small, with A < 2 ■ 1CT 3 . This fact is supported also 
by direct inspection of the largest pJW/, whose absolute values are three to four orders 
of magnitude smaller than the leading diagonal terms. The origin of the reduction of 
F p h in Fig. [7^ therefore lies in a modification of P n i, mediated by the off-diagonal terms 
of the density matrix. 

In Fig. [7] we show the population of the phonon bath Pi obtained in R 1 , panel (c), and 
in 5", panel (d). The changes mentioned above can be clearly seen: in Ri, for increasing 
7 the phonon population gets broader. In S an inversion of the population for the levels 
I = and / = 1 is present, with an even stronger sub-Poissonian F p h- 

The preliminary results show that there are regimes where dynamical coherent 
effects are relevant for the single phonon distribution. Much in the same way as for U, 
phonon distributions can be dramatically shaped by tuning 7. Strictly speaking, the 
results for 7 > 0.1 and ksT = 0.1u imply T L > k B T and therefore fall beyond the 
limits of the sequential tunneling regime for the calculation presented above. In this 
case we can expect that the additional coherent dot dynamics could lead to even larger 
deviations with respect to the RWA solutions. 

5. Conclusions 

In this paper we have analyzed the out of equilibrium properties of an unconventional 
single mode phonon bath, coupled to a quantum dot, in terms of its phonon distribution 
and the respective Fano factor. We have described the combined evolution of the dot 
and the bath degrees of freedom by means of a GME. In the regime of validity of the 
RWA, corresponding to a fast bath motion, we have shown that the out of equilibrium 
phonon bath markedly differs from a thermal one. The transport properties of the dot 
in the presence of a thermal or an unequilibrated single bath mode are qualitatively 
different. In the case of a single occupancy of the quantum dot, we have identified 
the transport regimes and the range of asymmetries A and dot-bath coupling strengths 
A for which a sub-Poissonian phonon distribution appears. We have interpreted such 
distribution in terms of a peculiar, selective population of the bath modes. Double 
occupancy of the dot may act as a filter on the phonon bath distribution, eventually 
turning F p h < 1 in regions where a super-Poissonian distribution would be attained in 
the case of a single dot occupancy. The crossover from an unequilibrated phonon bath to 
a thermal regime has been investigated by coupling the single bath mode to an external 
dissipative environment. For increasing damping strength, the sub-Poissonian phonon 
bath is destroyed and the latter tends to a thermal one. The crossover from a sub- to 
super-Poissonian bath has been analyzed as a function of A and A. We have found that 
F p h < 1 survives up to relatively strong damping when A ph 1 and A ~ 1. Finally, we 
have shown that in the case of uoo not much larger than the tunneling rates, the RWA 
is no longer valid and the solutions of the full GME display strong modifications when 
the coherent dynamics of the bath is fully taken into account. In particular, we have 
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shown that although off-diagonal density matrix elements remain small in the coherent 
regime, their coupling to the diagonal elements gives rise to significant modifications of 
the latter. These results show that by suitably tuning the quantum dot properties, it 
is possible to obtain a variety of unconventional single phonon baths characterized by a 
wide range of phonon distributions. A detailed analysis of the coherent regime and the 
possible extension to out of equilibrium baths with two or few more phonon modes is 
deferred to future work. 
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